Universality and saturation of intermittency in passive scalar turbulence 
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The statistical properties of a scalar field advected by the non-intermittent Navier-Stokes flow 
arising from a two-dimensional inverse energy cascade are investigated. The universality properties 
of the scalar field are directly probed by comparing the results obtained with two different types 
of injection mechanisms. Scaling properties are shown to be universal, even though anisotropics 
injected at large scales persist down to the smallest scales and local isotropy is not fully restored. 
Scalar statistics is strongly intermittent and scaling exponents saturate to a constant for sufficiently 
high orders. This is observed also for the advection by a velocity field rapidly changing in time, 
pointing to the genericity of the phenomenon. The persistence of anisotropies and the saturation 
are both statistical signatures of the ramp-and-cliff structures observed in the scalar field. 



Ramp-and-cliff structures are a characteristic feature 
of fields, like dye concentration or temperature, obeying 
the passive scalar equation (see, e.g., Refs. 0J^]) : 

(1) 



d t T(r, t) + v{r, t) • VT(r, t) = nAT{r, i), 



i.e. advected by the velocity v and smeared out by 
the molecular diffusivity k. Scalar gradients tend in- 
deed to concentrate in sharp fronts separated by large 
regions of weak gradients (see Fig. 1). The experimen- 
tal evidence for ramps and cliffs is long-standing and 
massive Furthermore, numerical simulations in- 

dicate that scalar structures are not mere footprints of 
those in v and appear also for synthetic flow |7,||]. The 
presence of ramp-and-cliff structures raises some impor- 
tant issues about scalar turbulence and its intermittency 
properties. Following Kolmogorov's 1941 theory, it is in- 
deed usually assumed that turbulence restores universal- 
ity, i.e. independence of the large-scale injection mecha- 
nisms, and isotropy at small scales (see Ref. [||). The ev- 
idence for scalar turbulence is however that anisotropies 
find their way down to the small scales, manifesting in 
the scalar gradient skewness of O(l), independently of 
the Peclet number ||-^]. This is due to the preferen- 
tial alignment of ramp-and-cliff structures with large- 
scale scalar gradients, present in most experimental situa- 
tions. For structure functions S n (r) — ((T(r) — T(0)) n ), 
this persistence is revealed by normalized odd orders 
S% n +i/ S^ +1 decaying more slowly than the expected 
r 2 / 3 (see Ref. M). Is this experimentally observed behav- 
ior signalling that small scales are fully imprinted by the 
large scales and that the universality framework should 
be discarded altogether? This is the first issue, raised 
in Refs. jl]j5|], that we shall investigate in this Letter. 
The second is about the consequences of cliffs for high- 
order intermittency. Their strength candidates them for 
the dominant contributions to strong event statistics and 
the issue raised in Ref. 0] is whether structure function 



scaling exponents are then saturating to a constant for 
high orders n. 

Numerical simulations are an ideal tool to analyze the 
previous questions, allowing to probe universality, by 
comparing the results obtained with two different types 
of injection, and saturation, by gathering enough statis- 
tics to capture strong events. Here, we shall take for v a 
2D flow generated by a Navier-Stokes inverse energy cas- 
cade . Universality is then understood as dependence 
of scalar properties on the injection mechanisms for this 
fixed v statistics. The scalar is injected at large scales, 
comparable to those where the inverse cascade is stopped 
by friction effects, and its properties are investigated in 




FIG. 1. A snapshot of the scalar fluctuation field for the 
injection by a mean gradient. Colours are coded according to 
the deviation from the average value (from white to black). 
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FIG. 2. The 3rd and the 5th-order parallel structure func- 
tions for the injection by a mean gradient. In the inset, local 
scaling exponents dS n (r) / dlog r. The measured exponents 
are £3 = 1.25 ± 0.04, ( 5 = 1.38 ± 0.07, with error bars esti- 
mated from rms fluctuations of local scaling exponents. 

the energy inertial range (see Ref. |llj for details). There, 
the velocity is isotropic, scale-invariant with exponent 
1/3 (no intermittency corrections to Kolmogorov scaling 
jl2],[ll]]) and has dynamical correlation times (finite and 
free of synthetic flow pathologies discussed in Ref. 0). 

As for scalar injection, a first choice is naturally sug- 
gested by experiments, where it usually takes place via 
a large-scale gradient. We assume then, as in Refs. j^D, 
that the average (T) = g-r and we integrate the equation 
for the fluctuations 9 = T — g ■ r, i.e. ([j]) with a source 
term — v ■ g on the right hand side. A snapshot of the 
9 field is shown in Fig. 1. The presence of the gradient 
g breaks isotropy and allows for asymmetries and non- 
vanishing odd-order moments in the scalar statistics. The 
second choice is a more artificial random forcing f(r,t) 
added to ([l]). Its motivation is to produce an isotropic 
statistics, e.g. by taking / Gaussian, with zero average 
and correlation function (f(r, t) /(0, 0)) = 5(t)x(r / L). 
The scale L where the injection is concentrated is taken 
comparable to the velocity integral scale. The equations 
for the scalar are integrated in parallel to the 2D Navier- 
Stokes equation for about 100 eddy turn-over times by 
a standard pseudo-spectral code on a 2048 2 grid. In the 
runs presented in the following, the diffusive term is re- 
placed by a bi-Laplacian, but it was checked by another 
series of simulations that using a Laplacian gives consis- 
tent results, although on less extended scaling ranges. 

Let us first show that the persistence of anisotropies 
observed in experiments occurs also in our case. Odd- 
order structure functions vanish in the randomly forced 
case. In the shear case they do not, except for separations 
r _L g. For non-orthogonal r's, the scaling exponents do 
not depend on the direction r and in Fig. 2 we present the 
parallel structure functions, i.e. r aligned with g. The 
resulting third-order skewness S^/S^ 2 scales as r 25 , the 




5 r T/<(5 r T) 2 > 1/2 

FIG. 3. Pdf 's of scalar increments normalized by their stan- 
dard deviations for three separations r — 2.5 x 10~ 2 , 5x 10 -2 , 
10 _1 in the inertial range. 

2nd-order exponent being ~ 2/3 (see Fig. 4). As in the 
experiments, the skewness decay is slower than the ex- 
pected r 2 / 3 . Furthermore, here enough statistics is accu- 
mulated to give access also to the 5th-order. The persis- 
tence effect is now dramatic as S^/S^ 2 ~ r~ a - 2 increases 
at small scales. Intermittency generates of course an am- 
biguity in the normalization, e.g. S§/S^ 4 is decaying, 
albeit very slowly. This reflects the fact that scalar incre- 
ment pdf 's change shape with r and one should then be 
specific about which part of it is sampled and the choice 
of the observable representative of the anisotropy degree. 
It is however unambiguously clear that local isotropy is 
not fully restored at small scales and the quality of the 
scaling laws found here indicates that this is a genuine 
effect, not related to finite Peclet numbers. 

More insight into this breaking of full universality is 
gained by analyzing scalar increment pdf 's and moments 
of even order, which are non- vanishing for both types of 
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FIG. 4. Low-order even structure functions. Local scaling 
exponents are shown in the inset. The measured exponents 
are < 2 = 0.66 ± 0.03, C4 = 0.95 ± 0.04 and ( 6 = 1.11 ± 0.04. 



2 



10° 



10° 



10" 



10* 



10 u 



10"' 



shear 



2 
1.5 
1 

0.5 







^^^^ n= 










.Cn 


x x x i d 




n 



12 
10 



10 15 



10"' 



10" 



1 



FIG. 5. Scalar structure functions of orders 10, 12 and 14. 
The scaling exponents are shown in the inset. 

forcing. Fig. 3 shows that the pdf's for the two types of 
injections do not have the same shape (the same holding 
when symmetric parts are taken). In the shear case, the 
separations r have been taken along the diagonal direc- 
tions, at angles <fi = ir/A with respect to g. This choice 
is motivated by the application of the procedure devel- 
oped in Refs. JT^ JT^] to the 2D case and permits removal 
of the first subleading anisotropic contribution oc cos 2<f) 
to even-order moments. The fact that the pdf's have dif- 
ferent shapes implies that the adimensionalized constants 
C n in structure functions S n (r) — C n (er)"^ 3 (L/r) n ^ 3 ~^ n 
are not universal, as it was also explicitly checked by di- 
rect comparison. Conversely, in Fig. 4 it is shown that 
scaling exponents of even order moments are the same for 
the two types of forcing. For the pdf's this means that, 
although having different shapes, the curves are rescaling 
with r in the same way. 

The picture emerging from these results is as follows : 
structure function exponents Cn are universal, while con- 
stants, and thus the pdf's of scalar increments, are not. 
The difference between isotropic and anisotropic situa- 
tions is that the non-universal constants C^n+x in odd- 
order structure functions vanish by symmetry for the 
former case, while they generically do not for the lat- 
ter. Structure functions present anomalous scaling and 
there is no full restoration of isotropy while going toward 
small scales. This picture of universality is weaker than 
in Kolmogorov's 1941 theory, but coincides with the one 
emerged for intermittency in the Kraichnan passive scalar 
model [T^-p^] (see also Ref. Q]). The velocity used here 
has finite correlation times, scalar correlation functions 
do not obey closed equations, yet the universality proper- 
ties are the same. This points to a broader validity of the 
mechanisms identified for the Kraichnan model and it is 
likely that the same universality framework generically 
applies to passive scalar turbulence. 

Let us now discuss the consequences of cliffs for the 
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FIG. 6. On the left, the pdf's in Fig. 3 multiplied by r~ c °° . 
The upper curves are for the random forcing and the lower 
ones for the shear case. On the right, cumulated probabilities 
for scalar fluctuations to exceed \T rma vs r, with A = 2.5. 

intermittency at high orders. Their singularity strength 
suggests that the scaling exponents might saturate at 
large orders, i.e. Cn tends to a constant Coo for large 
enough n. Physical self-consistency for the survival of 
steepening strong fronts is demonstrated in Ref. p8[ , 
where saturation is shown to imply that dissipation pref- 
erentially spares the cliffs with the largest jumps. High- 
order structure functions in our simulations are shown in 
Fig. 5, together with the Cn vs n curve, compatible with 
saturation. The same holds for ratios of two moments vs 
r or one moment vs the other. Note that, for any finite- 
size field, there are orders where the moments start to 
be spoiled and some strongest single structure having a 
dissipation width will plausibly dominate the statistics, 
as in Burgers' equation |l9| . The convergence of the mo- 
ments was inspected by the usual test of checking that 
(S r T) u V decays before the pdf V{S r T) of the scalar 
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FIG. 7. As in Fig. 6, for different values of A in the ran- 
domly forced case. The curves compensated by r~^°°, with 
Coo = 1.4, are shown in the inset. 
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FIG. 8. Scaling exponents of the 2nd, 4th and 6th-order 
scalar structure functions for three different roughness expo- 
nents 7 in the Kraichnan model. 

increments S r T = T(r)—T(Q") becomes noisy. An alterna- 
tive observable more reliable than moments (as less sensi- 
tive to the extreme tails) is looking, for fixed 5 r T, at how 
V varies with the separation r. Saturation is equivalent 
to the pdf taking the form V{5 r T) = r c °° Q{5 r T/T rms ) for 
6 r T sufficiently larger than T rms = ((T - (T)) 2 ) 1 / 2 . The 
collapse of the curves r~'» °'P((5 r r) in Fig. 6 is therefore a 
signature of saturation and gives the unknown function 
Q. In Fig. 7, we plot the cumulated probabilities f£° T V 
vs r for various S r T and the parallelism of the curves is 
again the footprint of saturation. Explicit evidence for 
the universality of Coo is provided in Fig. 6. 

The physical origin of cliffs resides in the Lagrangian 
structure of (0), i.e. in the fact that particles are pas- 
sively transported by the velocity v. In regions where 
velocity gradients are sufficiently persistent in space and 
time, widely spaced particles tend to approach and gen- 
erate the observed abrupt variations of the scalar field. 
This suggests that, even though quantitative aspects, 
such as the order of saturation or the value Coo, de- 
pend on the choice of v, the saturation phenomenon itself 
should occur for a wide class of random velocity fields. 
The Kraichnan model |^] is unfavorable for saturation 
because of the short velocity correlation time. Despite 
this, for large dimensionalities of space, saturation ana- 
lytically follows from an instanton solution |2(J. For the 
3D case, saturation was phenomenologically suggested 
in Ref. |21| and inferred from an instantonic bound in 
Ref. f22|| . Direct numerical evidence is provided by our 
3D numerical simulations whose results are presented in 
Fig. 8. Scaling exponents have been measured using the 
Lagrangian method presented in Ref. |^3| and (2 — j)/2 
is the spatial Holder exponent of v, as in Ref. @. The 
order of the moments needed to observe saturation is ex- 
pected to diverge for 7^2, while for 7 — > the action 
of large-scale gradients should favor close approaches be- 
tween particles. The order is thus expected to reduce 



with 7 and for the smoothest velocity in Fig. 8 satura- 
tion is indeed occurring already at the 4th-order and thus 
becomes observable. This confirms the physical picture 
of saturation due to the cliffs formed in the scalar field 
and the genericity of the phenomenon for scalar turbu- 
lence intermittency. 
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